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INFLUENCE  OF  FINITE  ION  LARMOR  RADIUS  EFFECTS  ON  THE  ION 
RESONANCE  INSTABILITY  IN  A NONNEUTRAL  PLASMA  COLUMN 

Ronald  C.  Davidson* 

Division  of  Magnetic  Fusion  Energy 
Energy  Research  and  Development  Administration 
Washington,  D.  C.  20545 

and 

Faculty  of  Science 

Hiroshima  University,  Hiroshima,  Japan 
Hwan-sup  Uhm 

Department  of  Physics  and  Astronomy 
University  of  Maryland,  College  Park,  Md.  20742  USA 

This  paper  investigates  the  influence  of  finite  ion  Larmor  radius 

effects  on  the  ion  resonance  instability  in  a nonneutral  plasma  column 

aligned  parallel  to  a uniform  axial  magnetic  field  Bq/£z-  The  analysis  is 

carried  out  within  the  framework  of  a hybrid  Vlasov-fluid  model  in  which 

the  ions  are  described  by  the  Vlasov  equation  and  the  electrons  are 

described  as  a macroscopic,  cold  fluid.  Electrostatic  stability  properties 

are  calculated  for  the  case  in  which  the  equilibrium  electron  and  ion 

density  profiles  are  rectangular  and  the  ion  distribution  function  is 
0 

specified  by  f ^-const. x6(Hi-wiP()-T^)G(vz) . The  resulting  eigenvalue 
equation  for  the  perturbed  electrostatic  potential  ^(r)  is  solved 
exactly  to  give  a closed  algebraic  dispersion  relation  for  the  complex 
eigenfrequency  w.  This  dispersion  relation  is  solved  numerically, 
and  it  is  shown  that  the  growth  rate  of  the  ion  resonance  instability 
exhibits  a sensitive  dependence  on  plasma  parameters.  For  example, 
finite  ion  Larmor  radius  effects  can  have  a strong  stabilizing  Influence 
for  azimuthal  mode  numbers  i>2,  particularly  when  the  equilibrium  self 
electric  field  is  sufficiently  weak.  For  the  fundamental  mode  (P“l), 
however,  stability  properties  are  identical  to  those  calculated  from  a 

[ 

macroscopic  two-fluid  model,  and  the  growth  rate  is  unaffected  by  the  value 
of  ion  Larmor  radius. 
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1 . INTRODUCTION 


One  of  the  most  basic  instabilities  that  characterizes  a nonneutral 

plasma^-  with  both  ion  and  electron  components  is  the  ion  resonance 
2-5 

instability.  In  cylindrical  geometry  (Fig.  1),  the  ion  resonance 

2 

instability  can  be  described  as  a two-rotating-stream  instability 

in  which  the  relative  rotation  between  electrons  and  ions  is  produced 

by  the  equilibrium  self-electric  field  E^(r)j|r.  Previous  analyses  of 

2-4 

this  instability  have  been  based  on  macroscopic  cold-fluid  models 
in  which  the  ion  and  electron  motion  is  assumed  to  be  laminar. 

Although  this  is  a reasonable  approximation  when 


we  expect  significant  modifications  to  the  stability  behavior  when 

r,  .J-R  . (Here  r,  , is  the  characteristic  thermal  ion  Larmor  radius, 

LI  ^ p Li 

and  Rp  is  the  radius  of  the  plasma  column.)  In  this  paper,  we 
investigate  the  influence  of  finite  ion  Larmor  radius  effects  on  the 
ion  resonance  instability  in  a nonneutral  plasma  column  with  rectangular 

electron  and  ion  density  profiles  (Fig.  2).  The  analysis  is  carried 

• • 

out  within  the  framework  of  a hybrid  Vlasov-fluid  model.  The  electrons 

are  described  as  a macroscopic,  cold  (T  -+0)  fluid  immersed  in  a 

uniform  axial  magnetic  field  B^e^.  On  the  other  hand,  to  allow  for 

the  possibility  of  large  ion  orbits  with  thermal  Larmor  radius 

comparable  to  the  radius  of  the  plasma  column,  we  adopt  a fully  kinetic 

model  for  the  ions  in  which  the  ions  are  described  by  the  Vlasov 

equation.  Such  hybrid  models  have  also  proven  quite  tractable  in 
6 7 8 

theta-pinch  ’ and  ion-layer  applications. 
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The  stability  studies  presented  here  assume  electrostatic  per- 
turbations with  infinitely  long  axial  wavelengths  (3/3z«0).  Moreover, 
the  analysis  is  carried  out  for  arbitrary  values  of  the  dimensionless 
parameter 


P 


For  a nonneutral  plasma  column  with  uniform  density,  it  is  also 

useful  to  introduce  the  dimensionless  parameter 

-1 

2ui 

5 = -fL  (1-f)  , 

U) 

ce 

which  is  a measure  of  the  characteristic  relative  strengths  of  the 

equilibrium  self-electric  force  and  the  magnetic  force  on  an  electron 

fluid  element.  Here  w is  the  electron  plasma  frequency,  id  is 

pe  r 7’  ce 

the  electron  cyclotron  frequency,  and  f=n^/n^  is  the  fractional 
charge  neutralization  provided  by  the  positive  ions.  In  the  present 
analysis,  the  parameter  5 and  the  fractional  charge  neutralization 
are  allowed  to  span  the  range  of  values 


and 


0<S<1  , 


0<f<l  , 


where  <$“1  corresponds  to  the  maximum  allowed  charge  density  for 
radial  confinement  of  the  equilibrium  configuration. 

The  outline  of  this  paper  is  the  following.  In  Sec.  2,  we 
describe  the  hybrid  Vlasov-fluid  model  (Sec.  2. A)  and  summarize 

the  equilibrium  formalism  (Sec.  2.B)  for  general  electron  density 

0 q 

profile  ng(r)  and  equilibrium  ion  distribution  function  of  the  form 

n n 


/ 


Here,  v^  is  the  axial  velocity,  is  the  perpendicular  energy 
(Eq.  (9)1,  P()  is  the  canonical  angular  momentum  [Eq.  (10)],  and 
const.  In  Sec.  2.C,  equilibrium  properties  are  calculated  for  the 
case  where  the  electron  and  ion  density  profiles  are  rectangular 
(Fig.  2)  and  the  ion  distribution  function  Is  specified  by 
[Eq.  (18)1 


0 fn0ml 

fi(*’*>  “ ~2 IT  MH1-.olP0-T1)G(v2) 


where  f,  n^  and  are  constants.  Electrostatic  stability  properties 
are  discussed  in  Secs.  3 and  4.  In  Sec.  3. A,  the  general  eigenvalue 


equation  [Eq.  (49))  Is  formulated  for  arbitrary  n (r)  and 

e 


VH.i  -co^P^,  v2) • In  circumstances  where  the  perturbed  ion  density 


corresponds  to  a surface-charge  perturbation  (at  r**lO  , the  eigenvalue 


equation  (49)  is  then  solved  in  Secs.  3.B  and  3.C  for  the  case  where 
,0 


f{  is  specified  by  Eq . (18)  and  the  electron  and  ion  density  profiles 
are  rectangular.  A striking  feature  of  this  analysis  is  the  fact  that 
the  required  orbit  integral  I [Eq.  (50)]  can  be  evaluated  in  closed  form 


[Eqs.  (59)  and  (62)]  for  general  values  of  the  parameters  r^/lt 


~2  2 

and  (2uipe/<nce)  (l-f) . Moreover,  the  resulting  eigenvalue  equation 


(61)  for  the  perturbed  electrostatic  potential  ^(r)  can  be  solved 


exactly  to  give  a closed  algebraic  dispersion  relation  ]Eq.  (66)]  for 
the  complex  eigenf requency  <o.  As  expected,  in  the  limit  whore  a-r^  ^-*0, 


Eq.  (66)  reduces  to  the  familiar  cold-fluid  dispersion  relation  ]Eq.  (70)] 

2,3 


previously  discussed  in  the  literature. 

The  general  dispersion  relation  (66)  is  an  algebraic  equation  of 


order  fc+3,  where  t is  the  azimuthal  mode  number.  In  See.  4,  a 
detailed  numerical  analysis  of  Eq.  (66)  is  presented  and  stability 
properties  are  Investigated  for  a broad  range  of  plasma  parameters. 
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It  is  found  that  the  growth  rate  of  the  ion  resonance  instability 

» -2  2 

exhibits  a very  sensitive  dependence  on  r,  ,/R  , w /<d  and  f. 

Li  p pe  ce 

For  example,  finite  ion  Larmor  radius  effects  can  have  a strong 

stabilizing  influence  for  mode  numbers  H>2  (see,  for  example. 

Figs.  10  and  14),  particularly  when  the  equilibrium  self-electric 

a 2 2 

field  is  weak  (w  /w  <<1  or  f close  to  unity) . For  the  fundamental 
pe  ce 

mode  (£»1),  however,  stability  properties  are  identical  to  those 

calculated  from  a macroscopic  two-fluid  model,  and  the  growth  rate  is 

unaffected  by  the  value  of  r,/R  . 

Li  p 

As  a final  note,  we  emphasize  that  the  present  stability  studies 
are  based  on  a fully  nonlocal  analysis  of  the  linearized  ion  Vlasov 
equation.  The  fact  that  the  dispersion  relation  can  be  obtained 
in  closed  form  is  a somewhat  fortuitous  manifestation  of  the  sharp- 
boundary feature  of  the  equilibrium  configuration  (Fig.  2)  and  the 
choice  of  equilibrium  ion  distribution  function  in  Eq.  (18). 


: 

5 

' 

\- 
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2 . THEORETICAL  MODEL  AND  EQUILIBRIUM  PROPERTIES 


2. A Theoretical  Model 


In  the  present  analysis,  the  electrons  are  treated  as  a cold 


(T  -*-0)  fluid  immersed  in  a uniform  axial  magnetic  field  B„e  . 
e & CKz 


Within  the  context  of  the  electrostatic  approximation  (J^-B^e^  and 


VxjE^O) , the  equation  of  motion  and  the  continuity  equation  for  the 


electron  fluid  can  be  expressed  as 


(~p  + V • V)V 
3t  ^e  ^e 


e 

m 


(-v* 


V xB„e 
"ve  O'vz 


) • 


(1) 


tr  n + V • (n  V ) =0  , 
3 1 e e've 


(2) 


where  £(£,  t)=-V<f>(£,  t)  is  the  electric  field,  n^(x,t)  is  the  electron 
density,  V^C^.t)  is  the  mean  electron  velocity,  and  -e  and  m 
are  the  electron  charge  and  mass,  respectively.  In  Eq . (1),  the 


spatial  variation  in  Bq  is  neglected  (low-beta  approximation) 


To  allow  for  the  possibility  of  large  ion  orbits  with  thermal 

Larmor  radius  comparable  with  the  radius  of  the  plasma  column,  we 

•*** 

adopt  a fully  kinetic  model  in  which  the  Ion  distribution  function 
f.(£,v,t)  evolves  according  to  the  Vlasov  equation 


h , 8 e / „ ^*Bo£z\  3 \ . 

[ i 


t)  - o , 


(3) 


where  +e  and  m.  are  the  ion  charge  and  mass,  respectively.  In  Eqs. 


(1)  — (3) , the  electrostatic  potential  <J>(x,t)  is  determined  self-consistently 


from  Poisson's  equation 


= -4ne( 


d vf . -n  ) 
i e 


(4) 
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Kqu.it  Ions  (l)  - |-tf  constitute*  a closed  description  o!  tin*  equilibrium 
and  electrostatic  stability  properties  v» t t be  nonneutral  plasma  column, 
and  torm  t be  theoretical  basis  tor  the  subsequent  analysis. 


.'.h  deneral  Kquillbrlum  Properties 


An  equilibrium  analysis  ot  Kqs.  (I  )-(**>  tor  general  steady- 

state  l,l/jt**0)  profiles  proceeds  in  the  following  manner . As 

illustrated  in  Klg.  I,  we, introduce  cylindrical  polar  coordinates 

l r , 0 , * ) with  /-axis  correspond ing  with  the  axis  ol  symmetry; 

r is  the  radial  distance  from  tin*  z axis,  and  0 Is  the  polm 

angle  in  a plane  perpendicular  t o the  z axis.  Kor  azimuthal  lv 

symmetric  electron  equilibria  i .>  /,10**0  and  ,>/,iz“lfl  charact  er  ized 

by  n ir)  and  V |xf-Vl  trie,,  it  is  st  raight  forward  to  show  t rom 
e io  " ed  \0 

Kq.  (J)  that  the  functional  form  ot  electron  density  profile 
0 

n lr)  can  be  spec  Hied  arbitrarily.  Moreover,  the  deviation  t rom 

equilibrium  charge  neutrality  produces  a radial  electric  tield 
0 ti 

1 fri"-,!^  ,'r  that  Int  luenci's  the  a.*  limit  ha  1 mot  ion  ot  the  election 

r 

tluld.  It  follows  t rom  Kq . ill  that  equilibrium  force  balance  in  the 
radial  direction  can  ho  expressed  as  -m  V^'trl  ’ r-  -of1  l r 1 — «»V^  |rlb,  c, 

t'  t't*  |*  0 

or  equ  l va  l out  l v 


+10.  ao  “ I' 

e ce  e r ce 


where  “eb  m c is  the  electron  cyclotron  trequenev,  v i 1 — V |tl  i 

ce  t>  o o oO 

is  the  angular  velocity  ot  an  electron  fluid  element,  and  is 

the  angular  , frequency  defined  by 

.0 


cK 

i 

i B 


r»0  ’r 


it*i 


8 


Wo  note  from  Eq.  (S')  that  there  are  two  allowed  equilibrium 

values  of  ud  lr).  Solving  Eq.  (5)  for  u)  (r)  gives 
e e 


where  the  upper  sign  (w  "u>  ) corresponds  to  a "last"  rotational 

O t* 

equilibrium,  and  the  lower  sign  corresponds  to  a "slow" 

rotational  equilibrium.  Evident Iv,  in  the  limit  oi  weak  radial 

electric  field  with  |w,. '<<ui  . it  follows  from  Eq . 1 7 1 t hat  ,o 

r.  1 ' ce  e ce 

and  ut  -to,. . 
e r. 

Eor  the  ions,  anv  distribution  function  f'fx.vl  that  is  a function 
’ i v \ 

of  the  single-part ic le  constants  of  the  motion  in  the  equilibrium  tields 
Is  a solution  to  the  steady-state  id/dt«0)  ion  Vlasov  equation. 

For  present  purposes,  we  consider  the  elass  ot  rigid-rotor  ton  Vlasov 
equilibria  described  by 


.0 

t 

i 


! ■ V* 


v 1 


tSl 


where 


u)j*const . . 


v is  the  axial  velocity,  11,  is  the  perpendicular 


euergv 


m . 


irl 


IV  Is  the  canonical  annular  momentum 

s 

t’  - m,  r (v  fr,.'  , , thh 

t‘  i 0 el 

and  u>  *08  ,/w.c  is  the  Ion  cyclotron  frequence.  Once  the  functional 
Cl  'l 

0 

form  of  f v > is  pooit led*  the  equilibrium  electrostatic 

potential  ^\r^  can  bo  determined  so  \ f-oons  i st  ent  1 v from  the 

• 

st  eadv-st  at  o Poisson  equation 


r 


~~  r - 4iu* { jd  vf%.-w  P ,v  )-nl(r)]  , 1 1 1 1 

r ,>r  ,lr  I i 1 t t»  z e 

where  Jd  vt^-n*(rl  Is  the  equilibrium  ton  density. 

An  important  feature  of  equilibrium  distribution  functions  th.it 

depend  on  H,  and  1',  exclusively  through  the  linear  combination  H,-u>  P 
t*  ■‘•iO 

is  that  the  mean  a/lmuthal  mot  ion  corresponds  to  a rigid  rotation 
with  angular  velocity  w -const.  Del inlng  the  mean  a/lmuthal  velocity 

* . , 0 . f 1 0 , 1 0. 

of  the  ions  by  V^lri-ud  vv^t^Wi  d vt  t,  it  is  straightforward  to 
show  from  Eqs.  (8)-(,ll))  tlMt  (r)-u'  r for  the  class  ot  ion  equilibria 
described  by  Kq.  iSl.  In  the  equilibrium  and  stability  analysis  that 
follows,  it  is  useful  t o introduce  perpendicular  velocity  variables 
appropriate  to  the  rotating  frame  ot  the  ions.  Defining 


V - v . v , 
x xi 


Vv  " Vv"u’iX  • 


lor  equivalently  V -v  and  V -v  r>  . it  is  straightforward  to  show 
r r 0 0 l 

from  Eqs.  I and  Ol’l  that 


'"i  ..2 


Ui  - ->ir0  “ f + 


' > > 1 


where 


.>  V7-V~+V~«V~+V~,  and  i*(rl  is  defined  bv 
1 x \ rt’ 

o mi  : 

Vlrl  - ef'  (r) y (w  +u>  ui^lr 


Note  in  Eq . (id)  that  m^ is  the  perpendicular  kinetic  energy  in 
a frame  of  reference  rotating  with  angular  velocity  u> ^ , and  vlrl  is 
the  effective  potential  in  the  rotating  frame.  Substituting  Eq.  l L" 
into  Eq . 1^1,  the  equilibrium  ton  density  profile  n* .frl-  d vt*  can  be 

1 J 1 

expressed  as 


10 


n i ( r)  = 2 it 


dV,V, 


dv 


(T  '■  * tw'  \) 


(15) 


Moreover,  it  can  also  be  shown  that  the  equilibrium  pressure  tensor 
in  the  plane  perpendicular  to  is  isotropic  with  perpendicular 
pressure  P^t  (r)=nj  (r)!^  ( r)  given  by 


0 . . _0  , . „ 
ni(r)T  ( (r)  = 2ir 


mi  2 0 / mi  2 \ 

dViVi  | dv2  ~Y  Vifi  \ Y Vi  + vz' 


(lb) 

Equation  (lb),  in  effect,  determines  the  perpendicular  ion  temperature 
profile  (r)  in  terms  of  the  equilibrium  distribution  function  f^. 

2 . C Sharp-Boundary  Equilibrium 

The  formalism  outlined  in  Sec.  2.B  can  be  used  to  investigate 

equilibrium  properties  for  a broad  class  of  electron  density  profiles 

n^(r)  and  ion  distribution  functions  f?(H,-iD,P  ,v  ).  For  purposes  of 
e i 1 i 0 z 

analytic  simplification  in  the  stability  analysis  in  Sec.  3.B,  we 
specialize  to  the  case  of  a sharp-boundary  equilibrium  (Fig.  2)  in 
which  the  electrons  have  a rectangular  density  profile,  i.e.. 


°/  % 
n (r)  = 
e 


n^=cWnst . 


0<r<R 


R <r<R 
P c 


(17) 


In  Eq.  (17),  R denotes  the  radius  of  the  plasma  column,  and  r=R 
P c 

is  the  radial  location  of  a grounded  conducting  wall.  For  the  ion 
equilibrium,  we  make  that  particular  choice  of  f^  that  also  gives 
a rectangular  density  profile,  i.e.. 


0 ^nnmi 

li  ■ — wii-"1vVclV  • 


(IS) 


where  f and  T.  are  positive  constants,  and  (1(^1  Is  the  parallel 
velocity  distribution  with  normalization 


dv  G(v  )=1.  Substituting 
z z 


Eq.  (18)  into  Eq . (15),  it  is  straightforward  to  show  that  the  i( 


density  profile  can  be  expressed  as 


0,  ^ 

n1(r)  = 


fn^=const . 


i|)(r)<T 


*(r)>T. 


In  the  region  where  the  electron  and  ion  densities  are  constant  [Eqs. 
(17)  and  (19)],  tne  solution  to  the  equilibrium  Poisson  equation  (11) 
corresponds  to  the  parabolic  potential 


.0,  , i 2 

* (r)  = -n  w • 


where  <J>  (r=0)=0  has  been  assumed,  and 

2lTnnec  a.O 

“e  = (1_f)  = " d%-  = const-  (2^> 

is  the  ^°y^0  rotation  frequency.  Substituting  Eq . (20)  into  Eq.  (14), 
the  effective  potential  i|<(r)  can  be  expressed  within  the  plasma  column  as 


1 2 2 

'l'(r)  = j r , 


where  S2  is  defined  by 

2 2 + - 

fi  =to_  a)  . -w.-w.u ) =(io  .-io  . ) (u>,-u) . ) (23) 

E ci  l i ci  i i i i 

+ 

and  is  defined  in  Eq . (31).  We  note  from  Eqs.  (17),  (19)  and 
(22)  that  the  electron  and  ion  density  profiles  precisely  overlap  with 
a common  radi”«  R (Fig.  2)  provided  iJ)(K^)=f^,  i.e.  , provided 


in,  j 2 

T"rp  ’ 


which  relates  T and  the  radius  R of  the  plasma  column.  Making  use  of 
1 P 

Eqs.  (22)  and  (24),  the  ion  density  profile  in  Eq . (19)  becomes 


12 


| fn  , 0<r-R  , 

0 0 I1 

n1(f)  - I (2‘>) 

'0  , R <r<R  . 

P c 


Comparing  Eqs.  (17)  and  (2r>),  it  is  evident  that  n”(r)  = fn^(r) , where 
f«const. “fractional  charge  neutralization. 


We  make  use  of  Eqs.  (lb),  (18),  (22)  and  (24)  to  evaluate  the 
perpendicular  ion  temperature  profile  T()((r).  Some  straightforward 
algebra  gives 


TLL(r)  = Ti|1-r/Rp]  (26) 

for  0<r<R^.  Note  that  T^(r)  Is  a maximum  (T ^) for  r=(),  and  decreases 
to  zero  at  the  edge  of  the  plasma  column  (r=R^).  For  future  reference, 
it  is  useful  to  introduce  the  ion  diamagnetic  frequency  defined  by 

u’di”(c/eB0nir)(:'/3r)(n?Tii)'  Makln«  uso  of  Eqs.  (2S)  and  (2b),  we 
f ind 


u) 


di 


V 


eB 


o r~ 


const . 


(27) 


Substituting  Kqs.  (2  V)  and  (2  4)  into  Eq.  (2  7)  readily  gives 

2 

“i+“i“ci  = wci(Vudi)  • (28) 


i 


which  relates  the  (constant)  frequencies  (o^  and  Wjj.  Equation 

(28)  is  simply  a statement  of  radial  force  balance  (of  centrifugal, 
magnetic,  electric  and  pressure  gradient  forces)  on  an  ion  fluid 
element  for  the  choice  of  equilibrium  distribution  function  in  Eq.  (18). 
Equation  (28)  is  a useful  relation  since  it  indicates  that  the 
angular  rotation  velocity  w.  cannot  bo  specified  independently  of 
to  (!)„  and  io , . . For  example,  zero  mean  rotation  of  the  ions  (to  =0) 

l'  li*.  Cl  1 I 

is  consistent  only  if  w which  corresponds  to  an  exact  balance 


of  electric  and  pressure  gradient  forces  on  an  ion  fluid  element. 

On  the  other  hand,  for  complete  charge  neutralization  with  t-1  and 

to  “0,  It  follows  from  Eq.  (28)  that  toT+to.io  “to  ,<o,,-  II  (tor  example) 
E i i c 1 c i d l 

2 2 , „ 

the  inequality  to,  <<to  , is  also  sat  1st  led,  then  the  two  sol  tit  ions  lot 
d 1 cl 

Wj  can  he  approximated  by  anti  ui^-w  which  correspond  to 

"slow"  and  "fast"  rotational  equilibria,  respectively.  final lv,  for 
the  rectangular  density  profiles  prescribed  bv  Kqs.  (17)  and  (?'>), 
it  should  be  noted  from  Kqs.  (7)  and  (21)  that  the  mean  azimuthal 
motion  of  the  electron  fluid  corresponds  to  a rigid  rotation  with 

t 

iti  ^10  -const . 
e e 


2.l>  Ion  Trajector ies  for  Sharp-boundary  Equilibrium 

Of  particular  Interest  in  the  stability  analysis  in  Sec.  ).H 
are  the  ion  orbits  in  the  equilibrium  t leld  coni igurat ion  |see  Eq . (20) 


B°(x) 

r\-  '\ 


t ."») 


0 , 
K (x) 

v 


K°(r)e 
r 'f 


...  | 

i o T ( xe  +ve  ) 
«>  b > x v v 


In  Eq.  (2*1),  we  have  introduced  the  electrostatic  bounce  frequency  m 
defined  by 


2nn()e 


“b  “ Wl 


(1-0 


> •) 

It  electric  forces  dominate  magnetic  forces  (o~s'*<o‘  >,  then  the 

b c t 

perpendicular  ion  mot  ion  is  simple  harmonic  at  frequency  to  b*  1,1 
general,  however,  the  electric  and  magnet  lc  torces  may  be  comparable, 
and  the  subsequent  analysis  indicates  that  the  perpendicular  ion 
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motion  is  biharmonic  at  frequencies  10^  and  where 

2 


ci 


4w, 


U) 


ci 

2 


1 ± 1 + 


5)  )• 


ci 

1/2 


(31) 


1 ± 1 + 


The  ion  trajectories  required  in  the  stability  analysis  satisfy 


d*’ 


dt’  ’ 


(32) 


dv ' 

*\j  e 

dt ' m , 


„ . *’X,W 


«<*’> 


c ) 


where  x'(t,=t)=x  and  v'(t'=t)=v.  Substituting  Eq.  (29)  into  Eq . (32), 

f\j  % 'X/ 

we  f ind 


and 


z'  = z+v  (t'-t) 
z 


i-jx'(f)  ■ -V'<f)  +»ci^ ry'(t')  , 

at 

A‘ ^ 2d 


(33) 


(34) 


dt 


,2  y ' (t ' ) =-v’(f)  -cl^rx-(t-)  . 


Equation  (34)  can  be  solved  exactly  to  give  the  particle  orbits  in 
the  plane  perpendicular  to  Bg£z-  In  this  regard,  it  Is  useful  to  make 
use  of  the  velocity  variables  (V^.V^)  in  the  rotating  frame  defined  in 
Eq.  (12).  Moreover,  we  introduce  the  polar  velocity  variables  (Vx,4>) 
in  the  rotating  frame  defined  by 

v +u).y  = V = V. co si  , 
x y x 1 


(35) 


v -u).x  * V = V,sin*  . 

y J y 1 


j 
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In  solving  Eq.  (34),  we  also  make  use  of  the  fact  that  the  Cartesian 
coordinates  (x,y)  are  related  to  the  polar  coordinates  (r,0)  by 
x=rcos6  and  y=rsin0.  Integrating  Eq.  (34)  with  respecr  to  t',  and 
expressing  the  resulting  solutions  for  x'(t')  and  y'(t')  in  terms  of 
( VA , <t> ) and  (r,0),  we  find 


x'(t)  = — j (V,  [sin(^T)-sin(*+u  t)] 


V“i 


(36) 


+ r(a)^-u)^)cos(0+a)^T)-r(wi-o)^)cos(0+w^T)  } , 


and 


y’(T)  = - ^ {Vi[cos(4.+a).T)-cos(^T)  ] 


Vi 


(37) 


+ r(u)j-<D . )sin(0-hu+T)  - r(to^-w^)sin(0+u)iT) } , 


where  T=t'-t,  and  uj~  is  defined  in  Eq.  (31).  It  is  evident  from 
Eqs.  (36)  and  (37)  that  the  perpendicular  ion  motion  is  biharmonic 
at  frequencies  w*  and  Shown  in  Fig.  3 are  plots  of  and 


versus 


Vi  - 2 2 

• Note  that  id.  asymptotes  at  +wu-w_. /2  for  . . 


ci  ~i  b ci b'  ci 

In  this  limit,  the  ions  are  effectively  unmagnetized  and  exhibit  simple 

harmonic  motion  at  frequency  w,  in  the  electrostatic  potential  if>^(r)  = 

b 

(m1/2e)w^r2. 


1.  ELECTROSTATIC  STABILITY  PROPERTIES 


lb 


I, A Ctntnl  Elaenavlue  Equat  ton 


In  this  section,  wo  linearize  Eqs.  (l)-(4)  assuming  electrostatic 

perturbations  about  the  general  class  of  axlsvmmct r Ic  equilibria 

described  by  t 1 1’^ , v ) |Kq.  lb)]  and  arbitrary  electron 

density  profile  i/\r)  (see  Sec.  l.B).  As  indicated  In  the  introduction 
e * 

the  present  analysis  assumes  flute  perturbations  with  3/3z-0,  so  that 
all  perturbations  have  spatial  dependence  onlv  on  ^■(x.y),  or  equiv- 
alently jtj_*(r,0).  In  the  elect  rostat  ic  approximation,  the  perturbed 
electric  field  is  AJ£(  Jtj. » 1 ) VAcS^>  ( x^ , t ) , and  Eqs.  (l)-(4)  can  he 
linearized  to  give 


(’  + u>  ~)iSV  - 1 -to  +2ui  )AV  . - v fit  , 

3t  e 30  er  ce  e e0  m 3r 

c 


( 18) 


(~  + Ut  i~)AV  , + 
3t  e 30  ert 


13  .9  , 

0 + , ( r ui  ) 

ce  r 3r  e 


e 1 3 

fiV  - -rA^.l  W) 

or  m r 30 


( + id  V ) An  + ; ~ (rn'vW  ) + AV  ■ 0 , (40) 

3t  e 30  e t 3r  e or  r 30  eO 


I 3 3 

— + v ■ — 

3t  ' 3 


+ 2 


♦ VV*1  . v UaM.o 


axr  1 


(4  1) 


m‘‘  IV*(?rt)l  * 3V  fi(^-Y.)  ’ 


ami 


— r + L ’-rrJiS*  “ -4ne(  d 'vAf  -An  ) • 

\r  ;,r  'r  r2  i 1 0 


(4.’) 


In  Eqs.  (l8)-(42),  AV  (j^.t)  1 s the  perturbed  electron  fluid  velocity, 
is  the  perturbed  electron  density,  l<i  1 j Ufi ' ' 1 ^ *s  • perturbed  ion 
distribution  function,  Jj'^(^)"-(34','/3r)y  Is  the  eqtti  I ibrltim  radial 


An  ^(^.t) 


electric  field,  (r)  satisfies  the  steady-state  Poisson  equation  ill), 
and  u>^(r)  is  the  equilibrium  angular  velocity  defined  in  Kq.  (7). 

To  simplify  the  right-hand  side  of  Kq.  (41),  use  is  made  of 

mfYj.  I*'1!'  (*01.  3^-ntj  ry((  [Eq.  (10)),  where  is  a unit  vector  in 

the  6-direction,  and  3f*?/ 3P .“-ui.  3f?/ 3H.  [Kq.  (8)).  The  linearized 

l o 1 1 1 

Vlasov  equation  (411  becomes 

{ If  m + ~ ($°<k>  + — r^)  * \ 

a o ' <43> 

’ e<.^"lJir^)  * VM^.t)  — fi(H*"uliW  • 


We  substitute 


, t ) “ ^i^)exp(-iu-t),  Imui'Q  , 


in  the  right-hand  side  of  Kq.  (4  *)  and  integrate  from  t'“-“>  to  t * = t 
using  the  method  of  charaet erist  les.  Neglecting  initial  perturbations 
and  noting  that  v and  df^/dl^  are  constant  (independent  of  t ’ ) 
along  a particle  trajectory  in  the  equilibrium  field  conf lgur.it  ion , 


we  t ind 


H<  r 

'S1  (^.v.t)  - e - dt  'expl-imt  ' ) 


•)!-»,  Jjt  .'(leu'll}  . 


where  ^j(t')  and  1,0  t,u'  perpendicular  particle  trajectories 

in  the  equilibrium  fields  and  /dr)j^  |Eq.  (V)|, 

and  >f*  (r)  satisfies  the  general  equilibrium  Poisson  equation  (II). 
Making  use  of  v_[  • V^,S^(  ^j)  - (d/dt  1 ) >S$(  ^j)  to  integrate  by  parts  with 
respect  to  t‘,  and  changing  variables  to  i“t'-t,  Kq . (44)  can  also 


be  expressed  as 


where 
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3fl) 

v;  Tv7  1'^^ 


(4s) 


+ j diexp(-iu)t)  ^tw  - uij  ~r) 

In  writing  Kq.  (4S),  wc  have  made  use  of  Eq.  (11)  to  express 

fi(Hi‘a,iro,vz)  " fi  (t  va  + *<r)*  0 • 


( An ) 


and  hence  3fj/3HA  - (m.i^)  ^3f^/3VA,  where  m^V^/2  is  the  perpendicular 
kinetic  energy  in  the  rotating  frame,  and  i|>(r)  is  defined  in  Kq.  (la). 

We  now  assume  perturbations  ot  the  form  ) “if  ( r ) c xp  ( i (O') , 

l,'ne(i^)“n  j (r)exp(  i (0) , etc.,  where  ( is  the  azimuthal  harmonic  number. 
Some  straightforward  algebra  that  makes  vise  of  Eqs.  ( ISf-yiOl 
readily  gives 


Anon  , (rl 
e v 


1 3 

( 

V tL° 

3 

* I 

.2  to 

* — JE£  A 

r 3r 

2 

v V 

1* 

3i 

r **  v‘ 

t' 

> 

r :\(r) 

1 

- $£<-*  +2ui  1 

3r  1 - ee  e 

ui-  (ui 

Tr) 

(a  7) 


2 0 2 2 

where  m (r)  Ann  (r)e  /m  , v"(r1  is  defined  bv 
pe  e i'  e 


v (r)  ” (io  -2ni  ) 
e i'  e i' 


ui  — ( r ui  3 

ee  r 3r  e 


- (u'~  v'u'  3' 


t ASf 


and  ui(r)  is  determined  from  Eqs.  (7)  and  (11>.  Moreover,  substituting 
Kqs.  (AS)  and  (A7)  into  Kq . (42),  the  linearized  1'oisson  equal  ion  can  be 
expressed  as 

if/  “C..\  i . ) i - / “C..V 

■) 


1 3 
r 3r 


(|  ♦ V)  I;  iinj  - (, ♦ v)u, 

\ v r \ v 


cj(g)  i a 

r ..'-I'm  3r 


V (-m  +2u'  1 

2 ee  e 


4ne2  f 3 1 3fJ  f.  f°  - 

- — d v y- ♦(r)+i(u-tui) I dT$(r')exp [iH(6 '-0)-iwr  ] 

i ' v JL  i.  J — oo 

where  the  abbreviated  notation  4i(r)=^(r)  has  been  introduced  in  Eq. 

(49),  and  r'(i“0)»r  and  0'(t”O)=0.  Equation  (49)  is  the  eigenvalue 
equation  that  determines  $(r)  and  the  complex  eigenf requency  w for  the 
class  of  ion  equilibria  described  by  Eq.  (46)  and  general  electron 
density  profile  n^(r)  and  angular  velocity  profile  wg(r)  consistent 

with  Eqs . (7)  and  (11).  In  this  regard,  we  note  that  the  polar  velocity 
variables  (V^.ifi)  defined  in  Eq.  (35)  are  natural  variables  for  the 
integrand  in  Eq . (49).  Therefore,  in  the  subsequent  analysis  of  Eq . 


(49),  we  make  use  of 


> 2 it  ra 

dv  = d<J>  dVxV± 

J0  Jq  " 


dv 

z 


Moreover,  it  is  convenient  to  introduce  the  phase-averaged  (over  4>) 
orbit  integral  I defined  by 


I = i 


dx4>(r ')exp  [iH(9,-0)-iu)x]  . 


3 . B Eigenvalue  Equation  for  Sharp-Boundary  Equilibrium 


We  now  specialize  to  the  case  where  the  equilibrium  electron 
and  ion  density  profiles  are  rectangular  (Fig.  2) , as  discussed  in 
Secs.  2.C  and  2.D.  The  corresponding  equilibrium  ion  distribution 
function  is  given  by  Eq.  (18),  which  can  be  expressed  in  the  equivalent 


f°  - fTrt|rv'-ii(1-^))G<'’z)  • 

P 

where  use  has  been  made  of  Eqs.  (13),  (14),  (22)  and  (24). 
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Moreover,  making  use  of  Kq.  (SO)  and  3n  ( r) /Dr--n„iS  ( r-R  ) |Eq.  (17)1 

e U P 

the  eigenvalue  equation  (49)  van  he  expressed  as 


f 7r  ''  ('  f ) h - 7 ('  + r) «' > 


> nr)  ce 

r <o-tu> 


“ u» 

i>  r >h.-r  > 

- V u)  2 p 

e v 

o 


> •> 


Jv1v1 


Jf0 

livz  v,  .vvf  r^+Cui-euij > i 


where 


v (r)-(to  -2 io  ) -(..i-fui  ) ““const . , 

e ci'  e o 


io  (r) 
pe 


4lTnoe 

io  - = const.,  (Kr'-R  , 

pe  m p 


K <i'-K 
P ‘ 


Kvident  ly,  the  angular  frequencies  w(,  | Kq . (21)1,  l Kq . (10)1, 

[Eq.  (31)  1 and  ui  [Kq.  (7)]  are  constant  (independent  of  r)  for 
the  case  where  the  equilibrium  electron  and  ion  density  profiles 
are  rectangular  (Fig.  2). 

It  is  evident  that  the  perturbed  electron  contribution  to  the 

right-hand  side  of  Eq.  (52)  (the  term  proportional  to  3n(  / 3r=-n  ,<S(r-R  ) 

e 0 p 

in  Eq . (45))  is  equal  to  zero  except  at  the  surface  of  the  plasma 
column  (r»R^).  Moreover,  it  can  be  shown  that  Eq.  (52)  supports  a 
class  iif  solutions  in  which  the  perturbed  ion  density  (the  term 
proportional  to  jd  'v  V, ' ,)f^/dV,  . . . in  Eq.  ( 52)  | is  also  equal  to 
zero  except  at  r-R  . In  this  case,  it  follows  from  tin'  linearized 
Poisson  equation  (52)  that  the  electrostatic  potential  <j>(r)  has 


the  simple  form 
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$(r)  = Ar 


(55) 


inside  the  plasma  column  (O-KP^).  ln  Eq . (55),  A is  a constant, 
The  required  orbit  integral  ] [Eq.  (50)]  becomes 


iA 


f° 


J- 


diexpf-lC  9-iuiT  ] 


r2n 


— ^ [r ' exp(  10  *)]*’, 


2 71 


(56) 


where  we  have  interchanged  the  order  of  the  x and  <f>  integrations. 
From  Eqs.  (36)  and  (37) t we  express 


where 


[r'exp(ie')]*  = [x’(x)+iy'(x)]£ 


— ^a-[  ("0  V^exp  ( i<j>)+ra  (x ) exp  (i6)  ] ^ , 
(<ui”u>i) 


(^(x)  = exp  ( iw^x ) -exp  ( iu+x ) , 


a2(x)  = (a)1-w1)exp(iu)^x)-(u.-u)'|)exp(ia)”x)  , 

+ 

and  to1  is  defined  in  Eq . (31).  Substituting  Eqs.  (57)  and  (58) 
into  Eq.  (56),  it  is  readily  shown  that 

* x f0 

T - i — I dxexp(-iwx)  [(u>  -w“)exp(iw*x) 

(lO^-uO  ' -«  11  * 

+ - £ 

- (u^-a)^)exp(iw^x)  ] 


(57) 


(58) 


(59) 


where  I (r)=Ar  . An  important  feature  of  Eq.  (59)  is  that  the  orbit 

A O 

integral  I is  indepe ''dent  of  perpendicular  energy  m.V^/2.  This  is 
a consequence  of  the  {articularly  simple  form  of  <Kr)  within  the 
plasma  column  [Eq.  (551]. 


Finally,  after  some  straightforward  algebra  that  utilizes  Eq . 
(51),  .he  ion  \e.loclt  integral  required  in  Eq.  (52)  can  be  evaluated 


to  give 


r"  r . 3f,  k 

2.  ] OV.vJ  d%  r W1  ■ ‘ ^ — 7r  "lM 
J0  '-•*  • V-1  'V1  (2T.  /m,)  ' 1 

1 l 


iS  ( r-R  ) 


( 2T  . / in . ) 
i l 


Substituting  Eqs.  (r''*)  and  (bO)  into  Eq . (S2),  the  linearized  Poisson 
equation  for  the  sharp-boundary  equilibrium  discussed  in  Secs.  2.1' 


and  2 . l>  can  be  expr 


esst'd  as 


r Yr  [r  ('  + f ) lr  *(r)  ~ ^ ('  + f ) ^r) 


. ui  u>  -2(0  (I)" , 

= *(r)  £ -4‘-  -lv  --  - l’1  R I' . (uO  6 ( r-R  1 . 

r j u>—  iu'  ' j p v p 

v e v . 

e l 

' 2 ' ' > 2 A •>  •> 

where  v 2T  /m  , ,o~  4nfn  e7m,  , (0~  4nn  e7m  , o'  eB,./m  c, 
t t i pi  0 i pe  0 e ce  0 e 

♦ 

u>  =,o  is  defined  in  Kq.  (7)  with  u>„  (2nn,.cc/B  .)  ( l-f)  = 
e e E 0 0 

, 2 2 

(Wpe^llVe^''^  lEq.  (20  1,  v~  is  defined  in  Eq.  (‘>3),  and  E , f,o)  is 
defined  bv 


. (u>)  - 


i(u>— (ui  ) 0 

1 - — — d i exp  ( — iu)T ) [ (m . — to . ) exp ( iw  r ) 

((0  , -10  . ) — " 

1 t 


+.  — . . < 
- (lo^-io . ) exp(  itOji ) | 


(0  -(0 , \ t . . (0-  <10.  / (0  . -(0 

1 i \ y 1 ' 1 / 1 i 

- + I „ m!(f-m)!  „ - , + ( + 

,o.-,o.  ' m=0  io-i,o.  -m(,o,-.o,)  \ ,o  , -,o 

it  l i i ' l i 


J.C  Dispersion  Relat  ion  tor  Sharp-Rpundarj 
Kquilib  t;  i um 


The  right-hand  side  ot  the  eigenvalue  equation  (bO  is  equal  to 


zero  except  at  the  surface  ot  the  plasma  column  (r=R  ).  Moreover , 


2 1 


rhe  i i 'ii.ii’i't  ion  Mi  > ..it  isi  i e;-  l lit'  vacuum  Poisson  equation, 

-1 


r • u |ro', I-  T/r  ) tf*M> 
solu,  i>  to  K(| . it>  . . te  expressed  as 


xc  pt  at  r*R  . Therefore,  tlu 
I' 


4> . (r  )*Ar  , t)  •-  r < K , 

P 


(b  31 


inside  the  column,  and 


1 f 0 C 

Il-R‘  /r-  ] 

^ (i)  - Ar^  — V. TT"  . R < r 

(t-irc/R2'i  P 


(b*) 


in  the  vacuum  region  between  the  sui race  of  the  plasma  column  and 


the  conducting  wall.  Note  that  4>(r)  is  continuous  at  r=R  , and 

P 


that  4^(r“K_)=»0.  The  dispersion  relation  that  determines  the  complex 


oigenfrequer cy  w is  determined  by  multiplying  Eq.  (bll  by  r and 


integrating  from  R ( 1 — t > to  R (1+t)  with  t:  *-0  . This  gives 
P P + 


R 


P I d t 


R 


r-R 


/ UC  > \ 1 * i ' 

1 + -~£  I ~— 

' v~  -,r  Ir-K 


U)  (i) 


_o 


o e* 


e 


O *■> 

U'  . K 

T((«l 

v . 
t 


(.  fc.  S) 


❖ fr=R  1 . 
P 


Substi'uting  Kqs . (b3l  and  (b4'  intc  Eq.  (o5)  and  rearranging  terms 
gives  the  dispersion  relation 


1 


I - ( R /R  1 
P c 


2< 


P!1 


1 ■> 
ui'lT 


+ r,w  . 


Obi 


2(w-tio  ll(u'-tu)  )-(u'  -2u'  11 

e e c e e 


vhere  . ,(td  is  defined  in  Eq.  (t>2),  to  =io  is  defined  in  F.q . (71  and  vise 
v e e 


has  been  mide  of  Eq . (331. 

'Equation  ob)  cut  be  used  to  investigate  the  stability  properties 


f • partial  I \ neutralized  electron  beam  for  a broad  range  of  dimensionless 


/ 2 

/ It 


aaramot ers  m , f,  etc.  \n  import  ant  parameter  that  characterizes 

pe  ce 


the  strength  ot  ion  kinetic  effects  is  t lie  ratio  of  charact  or  ist  tc  ion 


i 


Larmor  radius  (r,  “v,/w  ) to  radius  of  the  plasma  column  (K  ). 

Li  i c i p 

Making  use  of  Eqs.  (23)  and  (24),  it  is  straightforward  to  show  that 


rLi  Vi 

- s - ■■  - 

•>  o 1 

r~  H~u.r 

p p t 


( u>  -U)  . ) ( U).  -0)  , ) 

it  i i 


i wci 


where  is  defined  in  Eq . (31).  Evidently,  the  cold-ion  limit 

2 2 

(where  the  ion  motion  is  laminar  with  r“  <<R~)  corresponds  to  ion 

l.i  p * 

rotat  'ii  velocities  that  satisfy  w.-w.  or  w,  m,  . A careful  examination 

i i i i 

of  the  expression  for  r,(w)  (Eq.  (t>2)  ] shows  that 

t -2  2 i '2 

I i ^ at  . 

lia  _ I ^ r , i*w)  = ^ - — ioR) 

i ' i ""  Vi  2(i,>-v\j  ) | (u- C.o  ) + (w  +2w  ) ] 

i i c 1 l 


where  use  has  been  made  of  Eq.  (67)  and  u>  -u>  .+2u*  [Eq.  (31)  J. 

i i c l i 

+ 

Similarly,  in  the  limit  where  v , we  obtain 


u)*\R“ 

lim  + r{(u>)  = 

, *w  2 i v . 2 ( u>-  i io , 

i i v i i 


u)  , 

Ei 

, )((lO-tu)  ) + (ld  . +2ld.) 

i i c i i 


where  w -w  =w  ,+2».+.  Therefore,  making  use  ot  Eqs . (6R)  and  (6l>), 
i l c l i 

we  find  that  the  dispersion  relation  (66)  in  the  limit  of  zero  ion 
'2  ^ 

Larmor  radius  (r,  ,/R'  n))  reduces  to 
Li  p 


1 — ( K /R  )' 
P C 


J»£ + E-1 

- 1. ) [(oi-iii  )“(u  -2m  ) 1 2(m-£m. ) [(m-£m.)+(w  . +2u . ) 


e c e e 


i c i t 


where  (Eq.  (31) (.  Equation  (70)  Is  identical  to  the  familiar 

2 i 

dispersion  relation*’  obtained  within  the  framework  ot  a cold  two- 
fluid  model.  Throughout  the  remainder  of  this  paper,  Eq.  (70)  is 
referred  to  as  the  "reference  dispersion  relation"  (RPR).  In  order 
to  assess  the  influence  ot  finite  ion  gyroradius  on  stability  behavior 


it  will  be  useful  to  compare  the  stability  information  obtained  from 
Eq.  (6b)  with  that  obtained  f rom  Eq.  (70). 
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- ■>  ■>  _ 

Moreover,  In  tin*  laminar  cold-fluid  limit  whore  r *”  /R*"»0 , m «*m 

I.i  p i 1 

and  , tho  reference  dispersion  relation  (70)  can  he  expressed 


, 2 
ui 

pc 


( 7 S) 


1-(R  /R  ) 2(m-C<o  ) | (<o—  C to  ) - (w  -2m  ) | 

p i'  i*  e cc?  o 

' 2 

f ~ r-1 — 

2(io-£.u>, ) | (ui-Cu).  ) + (u  +2(ii  ) ] 

l 1 c i i 

where  and  to  are  defined  in  Eqs.  (11)  and  (71),  respectively. 

For  the  fundamental  mode  (£=11,  we  reiterate  that  the  linear  dispersion 

relation  (71)  is  identical  to  tho  reference  dispersion  relation  (7b). 

Therefore,  the  Vlasov-fluid  stability  properties  for  £=l  are  identical 

to  those  calculated  from  a macroscopic  two-fluid  model,  and  the  growth 

rate  is  una! tooted  bv  the  value  ot  P /R  . 

U p 

Hot  ore  proceeding  with  a detailed  analysis  ot  t lie  dispersion 
relation  (71)  and  comparison  with  the  reference  dispersion  relation 
(7S),  it  '.s  useful  to  summarize  the  limit  at  ion  t’n  plasma  parameters 
necessary  lor  existence  ot  the  nonneutral  equl librium  oont  i gurat i on . 
From  F(|s.  (71)  and  (72),  we  require 


pc 


1 1 - n - i , 


and 


i 1 


m , Jui 

f ) rv 

m 2 

c ut 

v*  O 


(i-n 


(7b) 


(77) 


tor  radial  cont  tnomont  o I the  oqu  i 1 ib  r i um.  The  inequality  in  Fq . (7b) 
assures  that  the  (repulsive)  space-charge  force  on  an  electron  fluid 
element  is  weaker  than  the  m.,,.net  ii  restoring  force.  Moreover, 

Fq.  (77)  is  equivalent  to  the  requirement  that  the  pressure  gradient 
force  on  an  ion  fluid  element  (in  t he  outward  direction,  since  )T^  /)r<0) 

In  Fig.  A , 


Is  weaker  than  the  confining  electric  and  magnetic  forces. 
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wo  make  uso  ot  Kqs.  (7b)  and  (77)  to  Illustrate  the  region  of  parameter 
'2  2 

space  (f ',m  /m  ) corresponding  to  existence  of  allowed  equilibria 

for  m./m  -183b  and  several  values  of  a-r,  ,/R  . Note  that  the 
i e 1.1  p 

equilibrium  electron  density  can  exceed  the  Brlllouln  flow  limit 

(u>*  /u'*-  ■O.S)  provided  there  is  sufficient  charge  neutralization  that 
pe  ce 

2 ' 2 

f Nl-0. 5(a)4"  ).  The  uppermost  curves  in  Fig.  4 are  calculated  from 

ce  pe  1 1 

the  inequality 


i (*  $ - .) 

m i V Rr,  / 


f V 1 - 


l_:  p 

2m*’  /m 
pe  ce 


which  Is  equivalent  to  F.q.  (77). 

The  growth  rate  and  real  oscillation  frequency  w "Rem  have 

been  obtained  numerically  from  l'q.  (7))  for  a broad  range  of  plasma 

% O •>  ^ 

parameters  f,  w*  /m*  and  r ./R  . We  now  summarize  the  essential 
pe  ce  1,1  p 

features  of  these  stability  studies.  The  analysis  is  restricted 
to  nonneutral  prot on-electron  plasmas  (ni^/m  -I83b),  and  the  growth 
rate  and  real  frequency  are  measured  In  units  of  t ho  lower  hybrid 
t requency 


..-(»'  , w 4) 

1,11  ce  cl 


Moreover,  throtighout  the  remainder  ot  this  section.it  is  assumed  that 


* - 0.1  . 


A detailed  study  o t the  Influence  ot  the  conducting  wall  locat Ion  on 


stability  behavior  is  summarized  in  Kef.  1 for  r,  . /R  "0. 

l.i  p 

1 2 

Stability  boundaries  lit  the  parameter  space  (f.m^  /<o‘  ) are 

pe  ce 

illustrated  In  Klgs.  '>-8  for  several  values  ot  v ^ ^ ranging  from  0 
to  2.1.  tn  Fig.  S,  the  solid  curves  correspond  to  the  st Ability 
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boundaries  (>=0)  obtained  from  Eq.  (73)  [or  Kq.  (75)]  for  r^/R^=0, 

and  several  values  of  azimuthal  harmonic  number  f.  For  a given 

'2,2 

value  of  < , the  region  of  (f,w  /u>“  ) parameter  space  above  the  curve 

pe  ce  

corresponds  to  instability  (y>0),  whereas  the  region  of  parameter 

space  below  the  curve  corresponds  to  stability  (3=0).  For  moderate 

■ 2 2 

plasma  densities  (w~  >0.01,  sav) , we  note  that  the  number  of 

pe  ce 

unstable  modes  increases  rapidly  as  the  fractional  charge  neutraliza- 
tion f is  increased  to  values  comparable  with  unity.  In  Fig.  6,  the 
solid  curves  correspond  to  the  stability  boundaries  obtained  from 
Eq . (73)  for  f|  j/R^=0.5,  and  C = 2,  A and  10.  It  is  important  to 
note  that  the  system  is  stable  in  the  region  of  parameter  space 
where  the  equilibrium  space-charge  field  Is  weak  [see  the  upper  left- 
hand  corner  of  Figs.  b(a)-6(c)).  Also  note  that  the  stable  region 
of  parameter  space  increases  in  area  as  the  azimuthal  mode  number  C 
is  Increased  [compare  Figs.  6(a)  and  6(c)].  In  Figs.  7 and  8, 

the  stability  boundaries  are  illustrated  for  r,  ,/K  =1  and  r.  ,/R  =2.5. 

l.i  p Li  p 

Evidently,  for  such  large  values  of  r,  ,/R  , the  region  of  (f,<o"  /to*-  ) 

Li  p pe  ce 

parameter  space  corresponding  to  allowed  equilibria  becomes  increasingly 
limited  by  the  constraint  in  Eq.  (78)  [see  also  Fig.  4],  However, 
the  qualitative  features  of  the  stability  boundaries,  and  their  dependence 
on  mode  number  C,  are  otherwise  similar  to  Figs.  5 and  6. 

The  dependence  of  stability  properties  on  fractional  charge 
neutralization  f is  further  illustrated  in  Figs.  9-11,  where  the 

normalized  growth  rate  and  oscillation  1 1 equency  K'r^u'i  q are 

„ ■>  2 ,92 

plotted  versus  f for  ui*"  /ui  =0.002  (Figs.  9 and  10)  and  m*-  /<o*‘  =0.5 

pe  ce  pe  ce 

(Fig.  11),  and  several  values  of  mode  number  ?.  We  also  assume 
fj  j/R^=0  in  Fig.  9,  and  r^  j/R^=0.5  in  Figs.  10  and  11.  In  Figs. 

9(b),  10(b)  and  11(b),  i^/ul  ^ is  plotted  only  for  the  ranges  of  f 
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corresponding  to  instability  (y^O).  Several  features  are  noteworthy 

in  Figs.  9 and  10.  First,  the  number  of  unstable  modes  increases 

rapidly  as  f is  increased.  Second,  for  8>2,  the  instability  growth 

rate  is  significantly  reduced  by  finite  ion  farmer  radius  effects, 

particularly  when  the  equilibrium  self  electric  field  is  sufficiently 

weak.  (For  example,  compare  Figs.  10(a)  and  9(a)  with  f approaching 

unity].  Third,  for  low  beam  densities,  we  note  that  the  real  frequency 

uir  exhibits  a nearly  linear  dependence  on  the  fractional  charge 

neutralization  f (Figs.  9(b)  and  10(b)]. 

The  dependence  of  stability  properties  on  beam  density  is  illus- 

■>  •) 

trated  in  Fig.  12  where  v/w...  and  w /ui, are  plotted  versus 

' l.H  r hli  pe  ce 

for  f-0.5,  r,,/R  -0.5  and  several  values  of  t.  As  shown  in  Fig.  12(a), 

Li  p 

the  absolute  maximum  growth  rate  (>,,  0.5b  w, ,,)  for  this  choice  of 

M l.H 

'2  2 

parameters  occurs  for  tm 2 and  ui*"  /ui  “0.38.  The  abscissa  in  Fig.  12 

pe  ce 

2 2 

extends  to  /w*"  *1.0  since  physically  allowed  equilibria  exist  for 
pe  ce 

“pe/wce-°,5/(l"°'5)“l  [Eq.  (7b)  and  Fig.  41. 

Of  considerable  interest  for  experimental  application  is  the 

*22 

stability  behavior  for  specified  values  of  f,  to*-  /m*  and  r /R  . 

pe  ce  Li  p 

Typical  results  are  summarized  in  Figs.  13-15,  where  and  (*>r/<*,lH 

* 2 2 

are  plotted  versus  mode  number  t?  for  /a)"  *0.01  (Figs.  I \ and  14), 

pe  ce 

Ci2  /ui2  -0.5  (Fig.  15),  r /R  -0  (Fig.  13),  r /R  -1.0  (Figs.  14  and  15), 
pece  Lip  Li  p 

and  several  values  of  fractional  charge  neutralization  f.  The  graphical 

results  are  presented  only  for  unstable  modes  with  >M).  For  sufficiently 

. ■> 

low  beam  densities  and  large  fractional  charge  neutralization  («*■  /ui*  - 

° pe  ce 

0.01  and  f-0.8  In  Figs.  13  and  14),  we  note  that  finite  ton  1. armor 


radius  effects  reduces  the  instability  growth  rate  and  the  number  of 


unstable  modes.  For  example,  for  r,  . /R  *0  and  f-0.8,  the  absolute 

l.i  p 

maximum  growth  rate  (y  -0.07  u ) occurs  for  f-12,  and  the  system  is 

M LH 


to 


unstable  for  4 in  the  range  l<t<18  [Fig.  13(a)).  On  the  other  hand, 
for  r ./R  =1.0  and  f=0.8,  the  absolute  maximum  growth  rate  (v  = 

LI  p M 

0.055  m^)  occurs  for  t= 9,  and  the  system  Is  unstable  for  1 < a < 1 A 
[Fig.  14(a)]  . For  high  beam  densities  (Fig.  15),  finite  ion  I.armor 
radius  effects  have  a negligible  influence  on  stability  behavior  (see 
also  Fig.  11),  although  there  is  a substantial  reduction  in  the 
number  of  unstable  modes  relative  to  the  low-density  case  (Figs.  13 
and  14) . 
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5.  CONCLUSIONS 

In  this  paper,  we  have  investigated  the  influence  of  finite  ion 
Larmor  radius  effects  on  the  ion  resonance  instability  in  a nonneutral 
plasma  column.  The  equilibrium  and  stability  analysis  (Secs.  2-4) 
was  carried  out  within  the  framework  of  a hybrid  Vlasov-fluid  model 
in  which  the  ions  are  described  by  the  Vlasov  equation  and  the  electrons 
are  described  as  a macroscopic,  cold  fluid  immersed  in  a uniform 
axial  magnetic  field  Moreover,  electrostatic  stability  properties 

were  calculated  for  the  case  in  which  the  equilibrium  electron  and  ton 
density  profiles  are  rectangular  (Ktg.  2)  and  the  equilibrium  ion 
distribution  function  is  specified  by  Kq.  (18).  In  circumstances 
where  the  perturbed  ion  density  corresponds  to  a surface-charge 
perturbation  (at  r*R  ) , a striking  feature  of  the  stability  analysis 
is  the  fact  that  the  required  orbit  integral  I [Kq.  (50)]  can  be 

evaluated  in  closed  form  |Kqs.  (5*1)  and  (b2)  ] for  general  values  of 

» '22 

the  parameters  r.  ,/R  and  (2io  /to'-  )(l-f).  In  addition,  the  resulting 
Li  p pe  ce 

eigenvalue  equation  (til)  can  be  solved  exactly  to  give  the  closed 

dispersion  relation  (bb).  A detailed  numerical  analysis  of  Kq.  (bb) 

was  presented  in  Sec.  4, and  it  was  shown  that  the  growth  rate  of  the 

loti  resonance  instability  exhibits  a very  sensitive  dependence  on 
* '22 

r,  ,/R  , »>*  /m"  and  f.  For  example,  finite  ion  Larmor  radius  effects 
Li  p pe  ce 

can  have  a strong  stabilizing  influence  for  mode  numbers  <^2  (see, 

for  example.  Figs.  10  and  14),  particularly  when  the  equilibrium 

'2  2 

self-electric  field  is  weak  (to  /to*"  <^1  or  f close  to  unitv).  For  the 

pe  ce 

fundamental  mode  (f«l),  however,  stability  properties  are  Identical 

to  those  calculated  from  a macroscopic  two-fluid  model,  and  the  growth 

rate  is  unaffected  bv  the  value  of  r,  ,/K  . 

Li  p 
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Fig.  10  Plots  of  (a)  growth  rate  y and  (b)  real  frequency 


versus  f [Eq.  (73)1  for  /u>^  “0.002,  m /m  =1836, 
M pe  ce  i e 

R /R  =0.5,  r, ,/R  =0.5,  and  several  values  of  i. 
pc  Li  p 


Fig.  11  Plots  of  (a)  growth  rate  y and  (b)  real  frequency 

* 2 2 

versus  f [Eq.  (73)]  for  ui  /to  *0.5,  m./m  =1836, 

^ pe  ce  i e 

R /R  *0.5,  r.  ,/R  =0.5,  and  several  values  of  t. 
pc  Li  p 

Fig.  12  Plots  of  (a)  growth  rate  y and  (b)  real  frequency 

versus  w"  /u>^  [Eq.  (73)]  for  f=0.5,  m /m  =1836, 
pe  ce  1 i e 

R /R  =0.5,  r,,/R  =0.5,  and  several  values  of  f. 
pc  Li  p 

Fig.  13  Plots  of  (a)  growth  rate  y and  (b)  real  frequency 

„ 2 2 

versus  4 [Eq.  (75)]  for  ui*"  /w  =0.01,  m /m  =183b, 

pe  ce  i e 

Rp/R  =0.5,  rj  an^  several  values  of  f. 

Fig.  14  Plots  of  (a)  growth  rate  y and  (b)  real  frequency 

versus  ? ]F.q.  (73)]  for  =0.01,  m /ra  =1836, 

M pe  ce  i e 

R /R  =0.5,  r, ,/R  =1.0,  and  several  values  of  f. 

pc  Lip 


Fig.  15  Plots  of  (a)  growth  rate  y and  (b)  real  frequency  w 

versus  4 [Eq.  (73)]  for  u>"  =0.5,  m /m  =1836,  R /R  1 

M pe  ce  i e p c 

0.5,  r,,/R  =1.0,  and  several  values  of  f. 
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